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ABSTRACT 

We have analyzed three XMM-Newton observations of the Seyfert 1 galaxy Mrk 509, with the goal to detect small variations in the 
ionized outflow properties. Such measurements are limited by the quality of the cross-calibration between RGS, the best instrument 
to characterize the spectrum, and EPIC-pn, the best instrument to characterize the variability. For all three observations we are able to 
improve the relative calibration of RGS and pn consistently to 4 %. In all observations we detect three different outflow components 
and, thanks to our accurate cross-calibration we are able to detect small differences in the ionization parameter and column density 
in the highest ionized component of the outflow. This constrains the location of this component of the outflow to within 0.5 pc of the 
central source. Our method for modeling the relative effective area is not restricted to just this source and can in principle be extended 
to other types of sources as well. 
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1. Introduction 

Active Galactic Nuclei (AGN) outflows are thought to play an 
important role in feedback processes, which con nect the growth 
of the black hole to the growth of the galaxy (Di Matteo et alJ 
I2005I) . However the strength of these outflows (kinetic luminos- 
ity) is unknown without knowing their precise distance from the 
central source. In the UV density sensitive line diagnostics have 
been used to constrain the density of th e outflowing gas and 
hence also the distance (Gab el et afll2003l) . Attempts have bee n 
made to do the same in the X-ray regime (Kaa stra et al.l l2004). 
but these results were affected by large uncertainties. Therefore 
in the X-ray band currently the best way to determine the dis- 
tance of the X-ray outflows is by using variability combined with 
high-resolution X-ray spectroscopy to consttain the density of 
the gas and then use this constraint to put an upper limit to the 
distance R. The main challenge with this approach is that it re- 
quires several crucial ingredients: 

- High-resolution X-ray spectroscopic observations to deter- 
mine the ionization state, column densities and outflow ve- 
locities of the warm absorber 

- Multiple observations of the same source 

- A source that varies on a suitable timescale and changes sig- 
nificantly in ionizing flux or in warm absorber properties 

If the source varies too fast, then it becomes difficult to obtain 
time-resolved spectra suitable to ionized outflow studies. If the 
changes in ionizing flux are small, then detecting changes in 
high resolution X-ray data is very difficult. There are several 
sources for which the variability approach has been successful, 
all with difficulties or with l arge changes in flux, however. The 
Chandra data of NGC 3783 dNetzer et alj|2003b [Krongold et alJ 
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I2005I) were spread out over several months (with distance lim- 
its between 0.2 and 25 pc), while the XMM-Newton observa- 
tions of this source (Be har et al.ll2003l) occurred during a grad- 
ual rising phase without clear flux extrema, with no changes 
in the absorber properties leading to a minimum distance of 
0.5 - 2 pc from the source, depending on the ionization com- 
ponent. In NGC 4051 t he flux dropped by a factor of five 
(ISteen brugge et al. 2009) leading to a 0.02 - 1 pc upper limit 
on the location of the absorbin g gas. For NGC 554 8 there was a 
3 year gap in the observations (Detmers et al. 2008) and the flux 
changed by a factor of five in the soft (0.2-2 keV) band, which 
resulted in an upper limit of 7 pc to the location of the warm ab- 
sorber. Also in NGC 3516 variability was used to constrain the 
location of the ion ized outflow (less than 0.2 pc from the source 
Netze r et al. 1 12001 . but the flux change was very large (50 at 1 
keV). 

In princi ple, the XMM-Newton EPIC-pn instrument 

dStruder et all 1200 ll) combined with th e XMM-Newton 
Refle ction Grating Spectrometer (RGS, Iden Herder et al.l 
2001) can be used to detect small changes in the column density 
of the ionized outflow. The RGS spectrum gives the detailed 
warm absorber structure such as outflow velocity, velocity 
width and ionization state. EPIC-pn with its larger effective area 
can detect small variations in the column density in response 
to small changes in the continuum flux. This requires that 
the relative calibration between EPIC-pn and RGS is better 
known than the observed variations in the column density and 
ionization parameter. 

To this end we have analyzed three different observations 
of the bright Seyfert 1 galaxy Mrk 509 (z = 0344 N H = 4.4 
x 10 24 nT 2 , Fishe r etalJll99 5: Murp hy et all 1 19961) . taken by 
XMM-Newton EPIC-pn and RGS, during different epochs in or- 
der to search for any variability in the outflow. We investigate 
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whether we can model the relative calibration between EPIC-pn 
and RGS in a consistent way, in order to take it into account 
when fitting the EPIC-pn and RGS spectra simultaneously. If 
successful, this approach allows for variability studies of weakly 
varying AGN, enabling us to constrain the outflow location also 
in these sources. Additionally this approach can be extended to 
other types of sources as well, as the method itself is source inde- 
pendent even though the exact model of the relative calibration 
will depend on the source type. We stress that it is also impor- 
tant, to take into account that the calibration of the instruments 
(e.g. effective area) has changed in time. 

Mrk 509 has been observed with high-resolutio n gratings be- 



fore, n amely in 2000 an d 2001 by XMM-Newton (Pounds et al 



2001 



2003). 



Paee et al.ll2003l) and Chandra HETGS (lYaqoob et al 



Smith et al. (2007) (hereafter S07) re-analyzed the RGS 



data and found evidence of three warm absorber components, 
as well as broad and narrow emission lines. We use this study 
as a starting point of our analysis of the thr ee new observations 
taken in 2005 and 2006. Using these data, ICappi et all d2009t) 
found evidence of a highly ionized and mildly relativistic out- 
flow. Clearly Mrk 509 is an excellent source to investigate in 
further detail. 

The most recent cali brations of the RGS and pn instruments 
(see lKaastra et al.l l2009, for the RGS) show that the uncertainty 
in the effective area of the RGS can be reduced to 3 %, while the 
pn absolute effective area calibration has an accuracy of 10 
Our aim is to obtain a more accurate relative calibration for RGS 
and pn. 

Sect. [2] describes the data reduction we have applied to all 
three observations. In Sect. [3] we describe the method we used 
to determine the cross-calibration differences between pn and 
RGS and how we model these. In Sect. l3.5l we present the results 
we have obtained for the cross-calibration. Sect. [4] contains the 
spectral analysis of all three spectra. We discuss our results in 
Sect. |5] and present our conclusions in Sect. [6] 



compare the results of the three individual observations with the 
stacked spectrum. 

The OM data were used to obtain fluxes in the UV band, 
which were used along with the pn data to construct the spectral 
energy distribution (SED) of Mrk 509. The fluxes in the UVW2, 
UVW1 and UVM2 filters were obtained by extracting the mag- 
nitudes from the PPS products of the OM data. These were then 
corrected for the extinction due to the Galactic absorption to- 
wards the source, using the relation between A^h and Ay from 
lPredehl&Schmitjq 995) and the relation between A,} and Ay, 
as given bv lCardelli et alj (1 19891) 

All the pn data were taken in small window mode with a thin 
filter. The background lightcurves were checked for flaring inter- 
vals and the good time intervals were determined to produce a 
filtered dataset. For the pn data we have first checked for possible 
pile-up effects. Using the epatplot task we analyzed the number 
of single and double events and found no signs of pile-up in our 
spectra. We also checked for X-ray loading, which can occur 
even when the countrate is below the pile-up limit. X-ray load- 
ing is the inclusion of X-ray events in the pn offset map which is 
produced at the start of the observation. This leads to additional 
and incorrect offset shifts for pixels where the X-ray loading oc- 
curs. All events in the subsequent observation associated with 
that pixel will be shifted to lower energies, depending on the 
energy of the incoming X-ray photons. A more detailed descrip- 
tion can be found in the EPIC calibration online documentation, 
document TN- 005(fl In order to counteract the effects of X-ray 
loading we excluded the central region from our analysis. We 
therefore took an annular extraction region centered around the 
source with an inner radius of 50 (2.5 arcsec) and an outer radius 
of 600 (30 arcsec) in detector coordinates. The background was 
extracted from an annulus of the same size outside the source re- 
gion and avoidi ng out of time event s. We use the SPEXEI spectral 
fitting package dKaastra et al.ll 19961) and we use C-statistics for 
fitting the spectrum. The errors are calculated for AC = 1. 



2. Data reduction 

The log of the three observations of Mrk 509 that we used is 
given in TableQ] We have used the SAS 9.0 software package to 
reduce all the data. We only use these three observations, as the 
S-N ratio of the first two RGS observations (SN of 3 and 4 in 
the 7-38 A respectively) is too low for this study due to a high 
background level. The S-N ratios of the three observations used 
here are 9, 7 and 8 respectively. We do not use the MOS data 
here, which are affected by pile-up and therefore their useful part 
has poorer S-N ratio as compared to the pn data. 

For the RGS spectra we use the rgspro$\ task with the source 
coordinates as obtained from SIMBAE0 (a = 311.04067 de- 
grees, 6 = -10.72353 degrees). After checking the background 
lightcurve for possible flaring intervals (observation 3 had one 
event with background level 0.25 cts/s and observations 5 had 
two events with background 0.4 cts/s, the mean net source count 
rate being on average 2.8 cts/s), we determine the good time in- 
tervals and rerun the rgsproc task now with the good time in- 
tervals selected (using a cutoff of 0.1 cts/s in the background 
lightcurve) and also with the keepcool keyword set to no, so cool 
pixels are also rejected. We also stacked all three observations 
with the rgscombine task for both RGS 1 and RGS2 in order to 



1 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 

2 http://xmm.vilspa.esa.eS/sas/8.0.0/doc/rgsproc/index.html 

3 http://simbad.u-strasbg.fr/simbad 



3. The relative calibration between pn and RGS 

3.1. Setup 

We only use data between 0.32 keV and 1.77 keV (this corre- 
sponds to the RGS band, i.e. 7-38 A). We follow a method in 
four steps, detailed below, to determine the cross-calibration be- 
tween pn and RGS. We first fit a model spectrum to the RGS 
data. All model parameters are frozen and the model is then ap- 
plied to the pn data. The ratios of predicted to observed pn count 
rates are then modeled to determine the relative calibration be- 
tween RGS and pn. The last step is to do a simultaneous fit of 
both the pn and RGS spectrum, taking the relative calibration, 
i.e. into account. As an additional check we created a fluxed RGS 
spectrum of each observation and follow the same procedure as 
above. 

3.2. Spectral Models 

We model all RGS spectra with a power-law and blackbody 
component, absorbed by three photoionized components that are 
modelled using the xabs model of SPEX. Also the cosmological 
redshift and Galactic absorption are taken into account. Broad 
emission lines are modeled with a Gaussian line profile, narrow 



4 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0050-l-0.ps.gz 

5 http://www.sron.nl/spex 
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Table 1. Mrk 509 Observations. 



Observation 


Date 


Obsid 


Duration 


Exposure RGS" 


Exposure pn° 


Flux" 










(ks) 


(ks) 


(ks) 


(W nr 2 s 


-') 


1 


18-10-2005 


0306090201 


86 


85 


58 


7.1 x 10- 


14 


2 


20-10-2005 


0306090301 


47 


47 


32 


7.2 x 10" 


14 


3 


24-04-2006 


0306090401 


70 


67 


43 


8.1 x 10" 


14 



" Elfective exposure time after filtering. 

b The flux is given for the pn 0.2-10 keV band. 



emission lines with a delta line profile. The flux at earth fx may 
then be written as: 



R(z)® 



fvh(Fph, Tpl) + /mbbCFmbb, ?mbb) + ^ Gj(Fi, At, cr ,) 



(=i 



+ Yjf A (A A ,F A ) X | |x ; (.V / .^.r ; .rr ; ; 



7=1 



A(N a ) 



We normalize the fluxes of the various components by their 
total flux in the RGS band. The power-law /pl is described by 
a normalization F pl and a slope Tpl, while the MBB /mbb is 
described by a normalization Fmbb and an effective temperature 
?mbb (see TableO. Each broad emission line G; is characterized 
by a normalization F,, central wavelength A,-, and width cr,, while 
the narrow line f A is characterized by a normalization F A and 
wavelength A A . Each absorption component X is characterized 
by a hydrogen column density Nj, ionization parameter out- 
flow velocity Vj and Gaussian velocity broadening cr,. The effect 
of cosmological redshift R depends on redshift z. The galactic 
absorption A depends on the hydrogen column density Nh- 

The widths cr,- of the two broad emission lines are kept fixed 
to the values determined from the optical observations and the 
central wavelengths Aj are fixed to the laboratory wavelengths 
(see Table 0). The high wavelength resolution of the RGS re- 
quires a detailed description of the wavelength dependence of 
the interstellar absorption. For this reason we use the SPEX 
model hot with a fixed low temperature kT — 5 x 10~ 4 keV and 
a fixed column density determined in earlier observations Nu = 
4.4 x 10 24 trr 2 (SectLQ. 

The ionization states of the photoion ized absorbers de- 
pend strongly on the assumed SED involved (Chakravortv et al. 
120091) . We have constructed the SED from the Mrk 509 observa- 
tions, using the pn and OM data to obtain the necessary flux 
points. The SED's of the three observations are almost iden- 
tical, so we use the SED obtained from the first observation. 
The photoionization balance calculations needed for our spec- 
tral modeling were based on this SED and performed using 
version C08.00 of C loud>B an earlier version is described by 
iFerland et ail d!998l) . The SED used is shown in Fig. Q] 

3.3. Method 1 

The first step is to obtain an accurate RGS fit for each obser- 
vation. This is done by fitting the spectral model described in 
the previous subsection (see Eq. [TJ to the RGS data. The second 
step is that we use this fit with all parameters fixed to predict 
the count rate Cm(PI) in the pn. We then compare these with the 
observed Co(PI) count rates in the pn by computing: 



t)2 




0.01 0.1 



10 100 1000 10' 10 J 
Energy (eV) 



Fig. 1. The SED used for the Mrk 509 observations. 
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The resulting values for R{PT) are plotted in Fig. IA.ll 

As the third step we calculate these ratios R(PI) using the 
SPEX model knak F. This is a piece-wise broken power-law, 
which consists of hinge points (A\, A%, .... ,Ak) and the corre- 
sponding transmission values (t\, £2, , fk)- We chose the hinge 

points by eye for the first observation, and keep these fixed for 
all other observations. The transmission values f,- are then fitted 
to the ratio R for each observation, leading to the values shown 
in Table [2] 

At this point we check the knak model F as follows: we com- 
pute F x /} and fold the result through the response matrix of the 
pn to predict the model count rate, and compare this with the ob- 
served pn count rate. Systematic deviations would indicate that 
the knak model must be revised, e.g. by redefining the A*. We 
find that the knak model given in Table |2]does not produce such 
systematic deviations (see Fig. IA.21 i. 

Now that we have a model F for the cross-calibration, we 
can fit the simultaneous RGS + pn spectrum. This is the fourth 
and final step. We create two sectors in SPEX, so that we can 
have two sets of models for the pn and RGS datasets, namely fx 
for the RGS and / pn = Ffx for the pn. We then couple the model 
parameters in such a way to take advantage of the strengths of 
each instrument. We let the RGS spectrum determine the high- 
resolution features, such as the fluxes F, and F A of the emission 
lines and the outflow velocity Vj and velocity broadening cr j of 
the ionized absorbers. This is because pn does not have the reso- 
lution necessary to accurately determine the detailed structure 
(broadening + outflow velocity) of the outflow. The pn spec- 
trum determines the parameters of the PL and MBB, and the 
column densities Nj and ionization parameters of the ionized 
absorbers (xabs). The resulting values are shown in Tables [3] H 
and [5] together with the values of the fixed parameters. 

We use this four step procedure for every observation, in- 
cluding the stacked one. 
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Table 2. Knak model parameters for the cross-calibration resid- 
uals. A value lower than 1 means the RGS flux is higher than the 
pn flux. 



A (A) 


1 




2 




3 


stacked 


7 


0.97±0.01 


1.00 


±0.01 


0.98 


±0.01 


0.99+0.01 


13 


0.98+0.01 


0.97 


±0.01 


1.01 


±0.01 


0.98+0.01 


15 


0.94+0.01 


0.97 


±0.01 


0.97 


±0.01 


0.96+0.01 


19 


0.96±0.01 


0.98 


±0.01 


0.98 


±0.01 


0.97+0.01 


23 


0.87±0.01 


0.92 


±0.01 


0.91 


±0.01 


0.90+0.01 


28 


0.92+0.01 


0.94 


±0.01 


0.91 


±0.01 


0.92+0.01 


34 


0.79±0.02 


0.78 


±0.02 


0.81 


±0.02 


0.80+0.01 


37 


0.90+0.03 


0.92 


±0.03 


0.94 


±0.03 


0.91+0.02 



10 20 30 

Wavelength (A) 

Fig. 2. Broad band RGS fit residuals for all observations. RGS 1 
data points are indicated by the black symbols, while RGS2 data 
points are indicated by the gray ones. The different observa- 
tions are marked by seperate symbols: triangle for observation 
1, square for observation 2, circle for observation 3 and cross for 
the stacked spectrum. 

3.4. Method 2 

As an independent check for each observation, we use the best- 
fit model fx for each observation to predict the RGS count rate 
Crgs and compute the model flux to be used for the pn f pn by 
multiplying fx with the ratio of the observed to predicted RGS 
count rates x fx- We then predict the pn count rate by 
folding f p „ through the corresponding pn response matrix and 
compare the result with the observed count rate. This proce- 
dure makes sure that any (time-dependent) astrophysical mod- 
eling problems are taken care of and the end result will be the 
pure cross-calibration residuals. We then compare the cross- 
calibration residuals we obtain with this method with the first 
method and whether the results are consistent with each other 
(see Fig. IA.ll ). 

3.5. Further checks 

To investigate the relative calibration further, we check for the 
presence of any weak deviations that might bias our RGS fit but 
would escape attention at full spectral resolution. Fig. [2] shows 
the residuals in the RGS band after rebinning the spectrum with a 
factor of 100. The typical residual fluctuations are at the level of 
3-4%. The third observation shows a discrepancy between the 
two RGS spectra. The RGS 1 data points are consistently below 
the RGS2 points, with a difference of about 2 %. This will have 
an effect on the total RGS flux, which we discuss later in more 
detail. 

Focusing on the pn spectrum, we load the best fit RGS model 
into the corresponding pn data set. If the relative calibration 
would be perfect, then the resulting fit quality of the pn spec- 
trum should be the same as for the RGS only fit. Instead we find 
something different. The resulting residuals are shown for each 
observation in Fig. IA.ll RGS and pn are in agreement with each 
other within 1-4% below 20 A. Above 20 A the RGS flux is on 
average 8-10% higher than the pn flux. All three observations 
and the stacked spectrum are in agreement with each other (apart 
from small fluctuations on the level of 2 %). 

As mentioned earlier we model these calibration offsets by 
a knak component. The best-fit parameters of this component 



are shown in Table [2] When applying this correction to each pn 
observation the fit improves dramatically. The resulting residuals 
are shown in Fig. IA.2l Some fluctuations on the order of 2 % can 
still be seen in the individual spectra, but our model adequately 
corrects for the large scale offsets. We chose the hinge points 
by eye and fix them for all other observations. The locations of 
these hinge points were chosen such, that they follow the trend 
visible in the pn residuals (see Fig. IA.ll in the Appendix). If we 
instead use a different number of points (e.g. 7 points, with fixed 
5 A spacing from 7 to 37 A), the residuals in Fig. IA. 21 change 
by no more than 1 -2 % We also performed additional tests with 
only 5 equidistant hinge points and the results are similar. 

Thus having obtained a model for the relative calibration for 
each observation, we perform a simultaneous fit of both the pn 
and RGS spectra, with the relative calibration differences taken 
into account by our knak model. 

4. Spectral Analysis 

4.1. Continuum 

We have modeled the continuum for all spectral fits with a 
power-law and a modified blackbody component. For all three 
observations and the stacked spectrum the continuum model pa- 
rameters are shown in Table[3] If we instead only use a powerlaw 
to model the continuum as was done in S07, the fits for all three 
observations are worse by AC = 75 per d.o.f. Using only the 
powerlaw as a continuum has a slight effect on the absorption 
components, as the powerlaw slope is softer in this case (r = 2.6 
instead of 2.3). However given the large change in AC, we use 
the continuum model of the powerlaw and modified blackbody 
from here on. The spectrum appears to be harder during the 3rd 
observation, while there is no significant difference between the 
first two observations (as expected, since they were taken only 
two days apart). The continuum parameters for the simultaneous 
fit are exactly the same as the RGS fits, but the errors are reduced 
due to the addition of the pn spectrum. 

4.2. Emission Lines 

Several broad and narrow emission lines are visible in the spec- 
trum (see Fig. [3] the O vn f narrow emission line (rest wavelength 
22.101 A lKeiivl[T987l) . as well as the broad O vn recombination 
line at 21.602 A (lEngstrom & Litzenl fl995l) and a broad Nvn 
emission line at 24.78 A dGarcia &Mack|[T96l . The oxygen 
lines were already detected by S07 (although they detect a broad 
intercombination line instead of a recombination line), but the 
N vn emission line is a newly detected feature. Table [4] shows 
the line parameters for the RGS only fit. The wavelength of the 
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Table 3. Continuum parameters for all observations. 



RGS 



Table 5. Warm Absorber parameters. 



RGS 



Observation 


F PL fl 


Tpl 


TMBB 


kT M BB (keV) 


1 


5.42±0.07 


2.30+0.03 


0.59+0.11 


0.092+0.006 


2 


5.73±0.08 


2.35±0.03 


0.67±0.17 


0.093+0.007 


3 


4.90±0.10 


2.06±0.04 


0.91±0.11 


0.129+0.006 


all 


5.44±0.05 


2.26±0.02 


0.68+0.09 


0.096+0.004 


RGS + pn 


Observation 


F PL fl 


Tpl 


TMBB 


kT MBB (keV) 


1 


5.41±0.05 


2.30±0.01 


0.61+0.06 


0.092+0.003 


2 


5.72±0.06 


2.35±0.02 


0.67+0.11 


0.093+0.005 


3 


4.94±0.06 


2.07±0.02 


0.87+0.05 


0.129+0.002 


all 


5.51±0.02 


2.27±0.01 


0.63+0.04 

r~- — TT^rnr 


0.095+0.002 



The unabsorbed flux in the 0.3 

Table 4. Emission line parameters for all observations, fluxes are 
corrected for Galactic absorption. 



Line 


Observation Wavelength 


Flux 


FWHM 






(A) 


(ph irT 2 s~') 


(A) 


O vn r 


1 


21.602 (f) 


3.0 


±0.4 


0.8(f) 


O vn r 


2 


21.602 (f) 


2.3 


±0.6 


0.8(f) 


O vn r 


3 


21.602 (f) 


2.3 


±0.5 


0.8(f) 


O vn r 


all 


21.602 (f) 


2.4 


±0.3 


0.8(f) 


N vn Lya 


1 


24.78 (f) 


1.0 


±0.3 


0.8(f) 


N vn Lya 


2 


24.78 (f) 


0.8 


± 0.4 


0.8(f) 


N vn Lya 


3 


24.78 (f) 


1.0 


± 0.4 


0.8(f) 


N vn Lycf 


all 


24.78 (f) 


1.1 


±0.2 


0.8(f) 


Ovnf 


1 


22.07 ± 0.01 


0.6 


±0.2 




Ovnf 


2 


22.07 ± 0.03 


0.2 


±0.2 




Ovnf 


3 


21.99 + 0.02 


1.0 


±0.5 




Ovnf 

1 tu. /-n : 


all 


21.99 + 0.02 0.30 


± 0.11 

1 * 





Comp 


Obs 


N H " 




V 








(10 24 m -) (10 -9 W m) (km s" 1 ) 


(km s~') 




1 


1.35+0.31 


0.60+0.18 


95+40 


- 230+70 


1 


2 


0.77±0.27 


0.34+0.20 


95 (f) 


- 230 (f) 




3 


0.61±0.21 


0.60+0.23 


95 (f) 


- 230 (f) 




all 


1.0±0.2 


0.6+0.1 


170+ 40 


- 120+ 60 




1 


12.0+2.1 


1.97+0.03 


80+20 


0+20 

-100 


2 


2 


11.4+2.5 


1.94+0.03 


80 (f) 


0(f) 




3 


6.1+1.7 


1.91+0.04 


80 (f) 


0(f) 




all 


10.5±2.0 


1.95+0.02 


70+ 10 


0+40 




1 


52+21 


3.27+0.09 


65+50 


-490+120 


3 


2 


26±12 


3.13+0.12 


65 (f) 


- 490 (f) 




3 


23±9 


3.06+0.09 


65 (f) 


- 490 (f) 




all 


SO™ 


3.20+0.08 


25+10 


- 290+40 


RGS + pn 




1 


1.36±0.23 


0.64+0.15 


95+40 


- 230+70 


1 


2 


0.79±0.22 


0.35+0.17 


95 (f) 


- 230 (f) 




3 


0.68±0.19 


0.63+0.20 


95 (f) 


- 230 (f) 




all 


1.2±0.2 


0.61+0.10 


170+40 


- 120+60 




1 


12.0+1.7 


2.01+0.02 


80+20 


"-no 


2 


2 


11.4+1.6 


1.96+0.03 


80 (f) 


0(f) 




3 


6.1+2.8 


1.92+0.04 


80 (f) 


0(f) 




all 


10.5+2.0 


2.01+0.03 


70+10 


0+40 




1 


51 + 14 


3.29+0.04 


65+50 


_ 475+60 

H/J -140 


3 


2 


25.4+6.2 


3.11+0.06 


65 (f) 


- 475 (f) 




3 


18.2+5.0 


3.01+0.06 


65 (f) 


- 475 (f) 




all 


69+so 
"*-m 


3.19+0.06 


25+10 


- 295+40 



The (f) indicates that the parameter was kept fixed. 
" Column density. 
* Ionization parameter. 
c r.m.s. velocity broadening 

d Outflow velocity, a negative velocity corresponds to a blueshift. 



O vn f line appears to shift between observations. However, this 
region suffers from bad pixels, which vary from observation to 
observation, so if one of the bad pixels falls on top of the loca- 
tion of the emission line, the fitted position shifts, depending on 
the surrounding pixels. Since these lines escape detection in the 
pn spectrum, we do not fit them in the simultaneous spectrum, 
but take the fitted values from the RGS only spectra. We have 
fixed the width of the broad emission lines to 0.8 A which corre- 
sponds to the width of the optical BLR (FWHM ~ 8000 km s _1 ). 
We also keep the wavelengths fixed to the rest-frame values. We 
do not detect any significant variability in the emission lines. 







i r 



i 1 



4.3. Warm absorber 

We detect three distinct warm absorber components, similar to 
what S07 have found by stacking two earlier RGS observa- 
tions of Mrk 509. We started out by using only one component, 
and then added components as the fit required. The second and 
third component are needed (A C = 50 and 40 respectively), but 
adding a fourth component does not improve the fit any further. 
The warm absorber in Mrk 509 consists of a low ~ 0.6), a 
medium ~ 2.0) and a high ~ 3.1) ionization component. 
The main absorption lines visible in the spectrum can be seen in 
Fig. [3] All warm absorber parameters are shown in Table|5] The 
values for the parameters are consistent for the RGS only and 
the simultaneous fit, but the errors on the latter are smaller due 
to the addition of the pn spectra. All parameters were set to a rea- 
sonable initial value (based on average values of warm absorber 



Fig. 4. The RXTE 2-10 keV lightcurve of Mrk 509, from 2005 
to 2007. Data were taken from the online RXTE mission-long 
data products. The diamonds indicate the exact date of the three 
XMM-Newton observations used in this paper. 



parameters found in the literature) and then freed to obtain the 
best fit. 



5. Discussion 

5. 1 . Warm absorber variability 

Comparing the warm absorber components for all three obser- 
vations we do not detect any significant variability in compo- 
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26 28 30 32 34 
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Fig. 3. The best fit RGS spectrum for the stacked observation with the absorption lines indicated. 
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nents 1 and 2. The high ionization component however shows a 
change in column density and ionization parameter between the 
first and the third observation. Changes in column density can 
be eit her due to motion a cross the line of sight or as in NGC 
5548 (Detmers et al. 200§1) due to recombination / re-ionization 
of the absorbing gas due to a change in the continuum flux. One 
has to be careful however in this case, as there is a correlation be- 
tween £ and column density, in the sense that a lower £ and lower 
column density can produce the same spectral fit as a higher 
larger column density case. In order to check for this effect we 
have performed also a fit where we have kept the column den- 
sity fixed to the column density of the first observation. All three 
observations have the same £ value for the high ionization com- 
ponent (log £ = 3.25). However for the third observations the fit 
in this case is worsened by AC = 37 for 1 d.o.f. with respect to 
the fit shown in Table[5] This indicates that the change in ioniza- 
tion parameter and column density between observation 1 and 
3 is significant. Using the RXTE monitoring data of Mrk 509 
taken with the PCA, we in principal can know what the con- 
tinuum flux history is of this source before and after the three 
XMM - Newton observations. However the gap during the first 
few months of 2006, prevent us from obtaining a more accurate 
timescale for the changes detected in the spectrum than the time 
between the 2nd and 3rd observation, namely half a year. If we 
make the assumption that the change in column density and ion- 
ization parameter we see are due to recombination of the com- 
ponent 3 gas due to the minimum, then we can get an estimate 
for the distance of this component using the following relation: 



nR 2 ' 



(3) 



£ is the ionization parameter, L the 1 - 1000 Rydberg luminosity, 
n the density of the gas and r the distance from the ionizing 
source. The 1 - 1000 Rydberg luminosity of Mrk 509 is 4.9 
x 10 37 W (as taken from the first observation). If we know the 
density of the gas then we can constrain the location of the gas, 
since £ and L are known from observations. In order to determine 
n, we must know the recom bination or ionization timescale of 
the ga s, which is given by (lKrolik&Krisslll995t iBottorff et al. 
2000): 



T KC (Xd = \aAXdn 



f(X M ) a r (X,_0 



f(Xi) a t (Xi) 



(4) 



where a r (Xi) is the recombination rate from ion X, + i to ion X t 
and f(Xi) is the fraction of element X in ionization state i. From 
the time between observation 2 and 3 we get an upper limit to r rec 
of 0.5 years. If we then take for example the FeXXII ion (since 
this is where the ionization level of iron peaks in column density 
for this component, see Fig|5]), we get a lower limit to the density 
n of 1.3 x 10 11 irT 3 . From this we can constrain the location of 
this outflow component to within 0.5 pc of the central source. 
This does not allow us however to distinguish between a scenario 
where the warm absorber is launched from the accretion disk 
(i.e. located at BLR distances) or from the dusty torus. Oxygen, 
neon, carbon and nitrogen are not produced in significant column 
densities by this high ionization component, so the only ions 
for which the column densities showed any significant change 
are Fe XXII to Fe XXV. If we perform the same excercise for 
the other ions, we obtain less strict upper limits (up to 4 pc for 
FeXXV). 

For the other two components which do not vary between the 
three observations, the non-variability can be due to two reasons: 




20 

Fe lonizotion level 



Fig. 5. Column densites for the different ionization levels of iron 
for the high (log £ = 3.2) component. The solid line represents 
the first observation, while the dotted line represents the third 
observation. 



The density is low, so that the recombination time of the gas 
is longer than the time between observations 3+4 and ob- 
servation 5, i.e. r rec > 0.5 years. 

The density is high enough to track the short time variability 
of the source, say on a few days timescale, as the fluxes of 
all observations are within 10 % of each other. 



5.2. Comparison with previous results 



Capp i et alJ ([2009) have found evidence of a highly ionized and 
variable outflow component using the same pn data we have 
used. In the RGS data there is no evidence of such an absorber 
with a mildly relativistic outflow velocity. As the main goal of 
this paper focusses on the cross-calibration of RGS and pn, we 
focussed our attention of the pn analysis only in the 7-38 A 
wavelength range. This highly ionized, fast outflowing absorber 
might represent an extreme part of the outflow described here or 
it might be disconnected given its high outflow velocity. 

If we compare the outflow components in our present spectra 
with those analyzed by S07, we see that the ionization param- 
eters are consistent with each other, although we find slightly 
lower values. The column density of the low ionized gas (log £ 
as 0.6) is a factor of 10 higher in the S07 observations. However 
the velocity broadening v they find is instead of the 105 km 
s -1 that we find. The main difference however is the outflow ve- 
locity. We find the highest outflow velocity (-425 km s -1 ) for 
the highly ionized (log £ ~ 3.1) component, while S07 find the 
highest outflow velocity for the component with log £ ~ 0.9. 

These differences may be due to a combination of several 
effects: 



L. 



Systematic uncertainty of 7 mA (50-200 km s ) in the RGS 
wavelength scale. 

A different continuum modeling in our paper compared to 
S07. 

Poorer S-N ratio in the S07 data. 

Sometimes an almost similar fit can be obtained by lowering 
the velocity broadening while increasing the column density, 
for components that are dominated by saturated lines. 
True time variability of the source. 
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Fig. 6. Comparison of the ratios R for both methods (see Fig 
IA.U for the stacked spectrum. The thin points are the R values 
for method 1 while the thick crosses are the R values for method 
2. 




10 20 30 

Wavelength (A) 



Fig. 7. Ratio of the pn residuals for each observation compared 
to the stacked spectrum. Observation 1 is indicated by the thin 
black line, observation 2 by the thin black line with dots and 
observation 3 by the thick black line. 



All this argues for a monitoring campaign where a high data 
quality will be obtained and the variations of the source are con- 
tinuously monitored. Such a campaign is now underway. 

5.3. Cross-calibration between RGS and pn 

This is a first attempt to model the relative calibration between 
pn and RGS in order to be able to detect small changes in out- 
flow properties. Our method for modeling the relative calibration 
is a relatively simple one, but a cursory glance at the differences 
between Figs. I A. 1 1 and IA.2I shows the clear improvements for 
the fit to the pn spectrum. If we compare the cross-calibration 
residuals for both the first method (RGS model for pn spectrum) 
and the second method (fluxed RGS spectrum folded through 
pn response), we see that apart from a difference for the first 
observation below 10 A, the residuals for both methods are al- 
most identical (see Fig. [6}. However there is a difference at long 
wavelengths, particularly from 22 A onwards. The same wave- 
like pattern which can be seen between 22 and 32 A in Fig [2] 
The difference is about 1.5% between the fluxed and the fitted 



data. The difference between the RGS fluxed and the fitted spec- 
trum can only be due to slight modeling inaccuracies, not due to 
any calibration offsets in the effective area of RGS or pn, as both 
methods use the same effective area. This means that by compar- 
ing the fluxed and the fitted residuals in the pn, we can separate 
the inaccuracies due to modeling and those due to uncertainties 
in the effective area calibration of RGS in the residuals seen in 
Fig. [2] These residuals are only apparent when the spectra are 
heavily rebinned. At higher resolution these residuals would not 
have been detected. From the comparison of the fluxed and fitted 
pn residuals we find that about 1 .5 % is due to modeling inac- 
curacies in the RGS fit. The strong wiggle in Fig. [6] between 22 
and 26 A with an amplitude of about 5 %, may be related to un- 
certainties in the pn model. Similar features can be seen in the 
residuals of for instance 3C273 orbit 1381 (see the calibration 
review tooQ). 

Both methods show the same general trend for the pn resid- 
uals. There is a sharp drop around 23 A, which is where the oxy- 
gen edge is located, possibly indicating that the modeling of this 
instrumental edge has some uncertainties. Comparing the resid- 
uals in Figs [7] and IA.2I between observations, we see that the 
general shape is exactly the same for each observation. There 
only appears to be a slight shift of ~ 3 % between observation 
1 and observation 3. The reason for this shift is the difference in 
flux in the RGS1 and RGS2 spectra of observation 3. Based on 
the comparison between the knak parameters, see Table [2j the 
RGS2 spectrum is the correct one. This means that in Fig. I A. 1 1 
the residuals in observation 3 should be shifted by 1 -2 % down- 
wards, which brings it into agreement with the previous obser- 
vations. The reason for this 2 % flux difference is unknown. We 
have checked for possible differences in exposure times or back- 
ground filtering criteria between the two RGS, and there is no 
difference between both RGS exposure times. The difference re- 
mains, independent of the exact criteria used for background fil- 
tering. As another check we varied the spectral extraction region 
for the third observation (90, 95 and 98 % of the cross-disperion 
PSF respectively). Also here the difference between RGS1 and 
RGS2 remains. At this moment we have no good explanation for 
this flux difference. It may however be a rare statistical variation 
at the 3.5 sigma level. 

If we compare the residuals with other cross-calibration in- 
vestigations, as shown on the calibration review tool website for 
XMM-Newton, we find identical features as in similar sources, 
like PKS 2155-304 and 3C273 for example, which show that be- 
low 0.5 keV the RGS flux is about 10 % higher than the pn flux. 
This gives us additional confidence that out method and our re- 
sults are correct and can be applied to other sources as well. 
However sources like Mrk 509 are the best to study the cross- 
calibration, as the pn data are not piled-up, but the source is still 
bright enough in the RGS to obtain a high-quality spectrum. 

The main application for this modeling of the relative cal- 
ibration is the increased sensitivity to changes in the outflow 
column densities and the ionization parameter the simultaneous 
fit provides. A comparison of Table [5] shows that performing 
a simultaneous fit improves the accuracy of the measurements 
(column density and ionization state) by almost a factor of two. 
Without correcting for the relative calibration, this would not 
have been possible. 

Our method has several assumptions however. The first is 
that the RGS data give us the best fit for the warm absorber and 
the emission lines. Also we assume that the RGS continuum fit 
gives the correct continuum. But since we are interested only in 

7 http://xmm2.esac.esa.int/cgi-bin/ept/preview.pl 
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the relative calibration, we do not draw any conclusions as to 
which of the two instruments is better calibrated in an absolute 
sense. We use the RGS to determine the best fit to the spectrum, 
because of the high-resolution needed for the precise determina- 
tion of the outflow properties. However as can be seen from Fig. 
[2] the RGS spectrum shows broad band fluctuations on the order 
of 3-4 %. As already mentioned earlier, this sets a limit on how 
good our relative calibration modeling can get. Also the exact 
knak model we have is at this moment strictly valid for sources 
with the same spectral type as Mrk 509 as there are indications 
that the relative calibration differences appear to be source and 
pn mode dependent This indicates that the redistribution of pn 
is partly responsible for the differences in the RGS and pn flux, 
since a different spectral slope will have an important effect on 
the exact redistribution of photons in the pn spectrum. This is 
also what has been reported to the XMM-Newton Users Grourfl 
By applying this method to other sources (including non-AGN 
sources) and looking at the exact differences between the ob- 
tained relative calibration residuals (i.e. the exact values of the 
knak model for example), gives us an opportunity to try to deter- 
mine the exact origin of the differences between RGS and pn. 

6. Conclusions 

We have investigated the cross-calibration differences between 
XMM pn and RGS. By modeling the cross-calibration differ- 
ences, which are consistent for all three observations and the 
stacked spectrum, we are able to improve the calibration for pn 
and RGS to within 4% of each other. With this improvement 
we can detect small variations in the outflow column density 
and ionization state more accurately. For Mrk 509 we have used 
this method in combination with the RXTE light curve to con- 
strain the location of the highly ionized component (log £ w 3.2) 
to within 4 pc of the central source. The other components do 
not show any significant variability, despite the source varying 
in flux between the 2005 and 2006 observations. With this im- 
provement to the relative calibration of the XMM instruments, 
we can use the advantage of simultaneous pn and RGS observa- 
tions to detect and model small changes in the outflow of AGN. 
This will allow us to expand outflow studies to sources which 
vary more slowly and with a smaller amplitude, thereby increas- 
ing the number of Seyfert 1 galaxies that can be studied in this 
way. This method can also be expanded to other types of sources, 
as the method itself is source independent. 
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Appendix A: Additional figures 
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Fig. A.2. The best fit residuals for the pn spectra after applying 
the knak model. From top to bottom: Observation 1, 2, 3 and the 
stacked spectrum. 
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Fig. A.l. The ratios R(PI) = c ° (p/) for all observations. The left-hand column shows the ratios for method 1, while the right-hand 
column shows the ratios for method 2. The locations of the knak hinge points are indicated by the triangles (The y-axis values for 
each point are not accurate, as the values shown in Table [2] are averaged between hinge points, the triangles shown here are just to 
guide the eye). From top to bottom: Observation 1, 2, 3 and the stacked spectrum. 



